Penguin and Box Diagrams in Unitary Gauge 

Xiao-Gang He* and Jusak Tandean''" 

Department of Physics and Center for Theoretical Sciences, 

G> ' National Taiwan University, Taipei 106, Taiwan 

O ' 



O 



o 

(N 



G. Valencia-'- 



I Department of Physics and Astronomy, 

Iowa State University, Ames, Iowa 50011, USA 



(Dated: September 20, 2009) 



Abstract 



We evaluate one-loop diagrams in the unitary gauge that contribute to flavor-changing neutral current 
j-^. (FCNC) transitions involving t-wo and four fermions. Specifically, we deal -with penguin and box diagrams 
arising -within the standard model (SM) and in nonrenormalizable extensions thereof -with anomalous cou- 
plings of the W boson to quarks. We sho-w explicitly in the SM the subtle cancelation among divergences 
(yQ from individual unitary-gauge contributions to some of the physical FCNC amplitudes and derive expres- 
' sions consistent -with those obtained using gauges in the literature. Some of our results can be used 
• more generally in certain models involving fermions and gauge bosons -which have interactions similar in 
, form to those we consider. 

o 
a^ 
o 

> 
X 



*Electronic address: hexgOphys . nt-u . ed-u . tw 
^Electronic address: jtaiideaii@yahoo.com 
■'■Electronic address: valencia@iastate.edu 



1 



In the presence of new physics affecting primarily the charged weak currents involving quarks, 
the W boson may have couplings to quarks beyond those in the standard model (SM). The effective 
Lagrangian for a general parametrization of the W boson interacting with an up-type quark U and 
a down-type quark D can be written as 

CuDW = -UY{9EPl + 9iPK)DW; + H.c. , (1) 

where g^^^ are complex coupling constants and Pl,r = i(l -F 75)- SM limit, these constants 

become 

9E = -^^UD, 5r = 0, (2) 

where g is the weak coupling constant and Vki are elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. 

Beyond the SM, the new couplings in Eq. (1) not only affect weak decays through tree-level 
interactions, but also modify flavor- changing neutral current (FCNC) transitions at one-loop level. 
In Ref. [1], we have studied several two- and four-fermion FCNC transitions induced at one-loop 
level by the new couplings in the charm sector and determined constraints on them from various 
processes. Here we provide the detailed derivation of the loop formulas summarized therein. Some 
of the results we present in this paper can also be applied to other models involving fermions and 
gauge bosons which have interactions similar in form to those we consider. 

In treating the loops, we adopt the unitary gauge, which is convenient in the absence of knowledge 
about the new degrees of freedom. Another advantage is that, since diagrams involving unphysical 
states are absent in unitary gauge, the number of diagrams to deal with is smaller than that in the 

gauges usually employed in the literature. Unlike in gauges, our calculations in unitary gauge 
produce divergences in some of the individual contributions. To treat them, we employ dimensional 
regularization with a fully anticommuting 75. We express the divergent part of the L>-dimensional 
integration in terms of the combination 

where 7^ is the Euler constant and ^ the mass scale that arises from the regularization procedure. 
For physical processes, the divergent parts should cancel in the total amplitudes in the SM limit 
and the results should agree with those obtained in R^ gauges. It is therefore instructive to see how 
the individual contributions in unitary gauge differ from those in the general R^ gauge. 

At the quark level, the loop-induced physical processes of interest in Ref. [1] are d — > d'7, 
d — > d'g, dd' — > vD, dd' — > and dd' — > dd' , where d and d' are down-type quarks, with 7, g, 

u, and i as usual denoting a photon, gluon, neutrino, and charged lepton, respectively. The relevant 
diagrams are displayed in Figs. 1 and 2, where the loops in unitary gauge contain only fermions 
and W bosons. The dd' vv and rfrf' — > t^i~ transitions also receive Z-penguin contributions 
and, for the latter, 7-penguin contributions as well. 

Some of our results are new, while some have been found before using unitary or R^ gauges. 
Specific aspects of our paper not previously available in the literature include unitary-gauge expres- 
sions for individual amplitudes corresponding to the separate diagrams in Fig. 1 contributing to 
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(a) (b) (c) 

FIG. 1: Diagrams contributing to amplitudes for d d'V*, with V being a neutral gauge boson. In all 
figures, straight lines denote fermions and the loops contain W bosons besides fermions. 
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FIG. 2: Box diagrams contributing to amplitudes for (a) dd' — > I'v or i'^i and (a,b) dd' — > dd' . 

d — i> d'V*, arising from the quark- 14^ interactions of Eq. (1) and SM-like couplings of the neutral 
gauge bosons (Z, 7, ^f) to fermions and the W . These particular results, presented in the Appendix, 
can thus be used more generally in models involving fermions and gauge bosons having interactions 
similar in form to those we discuss. Another new aspect of our work is the explicit demonstration, 
in unitary gauge, of cancelation between divergences from penguin and box diagrams contributing 
to dd' — vv and dd' —>■ i~^£~ . Moreover, perhaps more completely than earlier work in the liter- 
ature, we provide a set of unitary-gauge results for the FCNC transitions listed above induced by 
the interactions in Eq. (1), along with SM-like interactions of the neutral gauge bosons. We will 
compare our results with those in the literature, where available. 

In the following, we derive first the amplitudes for d d'V*, where d and d' are on-shell, but 
the neutral gauge-boson V is off-shell with momentum k that is small compared to m^r■ We have 
collected the resulting formulas corresponding to the diagrams in Fig. 1(a), (b), and (c) in the 
Appendix. Then, assuming that V has SM interactions with fermions, the couplings being given in 
Eq. (A7), we combine the three contributions. Thus we obtain for V = Z 



upon setting set A;^ 



9^1 
0, where 
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d 



(4) 



gtgf 



'±x z,x -r J- 
+ ■ — — Ina^q 



7a;„ A + —r^ ^ — ^ Inx. 



2(1-^,) 



(5) 



Zl 



gi gi zf^ i^q) + gi. gi. i^q) + gt gi ^q i^q) 
zl = gigtm^zf^i^q) + gigim^' z^{xq) + gtgim^M^q) 



(6) 
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13x - 2x^ + 7x- 



8x^ 



24(1 - x)3 



Zo{x) 



rSM, 



X] — 



12(1 

16x - 29x^ + 7x^ 
24(1 



2 x + x"^ - (4x^ - 6x^)c^ 
c,., H 77T 71 — Inx , 



4(1 

2a; - 



4(1 



Inx , 



5(4 + x + x^) - 4(20 - Six + 5x^)cl 5x - A{2x - ?,x^)cl 



24fl 



+ 



x] 



4(l-a;)3 



Inx 



(7) 



Evidently, the current (7^) terms are divergent, whereas the magnetic (o'^'^) terms are finite. Within 
the SM, the term containing zl has been calculated previously using unitary gauge in Refs. [2, 3], 
and our expression above agrees with the results found therein. 

In writing down M.'^^^d,'z*i have dropped finite terms (without A) that do not depend on the 
internal quark mass m^, or equivalently x^ — mym^. This is because such terms will be removed 
by the Glashow-Ihopoulos-Maiani (GIM) mechanism in the SM, after summing over q — u,c,t and 
imposing the unitarity relation V*^,V^^+V^,V^+V^*^,V^^ = 0, and by a GIM-like mechanism beyond 
the SM. For the same reason, we also drop m^-independent finite terms in the next two amplitudes 
for d — * d'V*, as well as in the box-diagram contributions below. 

For V = 7, after combining the three contributions (a,b,c), we find that the /c-independent terms 
cancel completely, but that terms of first and second orders in k do not, which we keep, namely 



ee, 



d' 



d 



(8) 



where 



d -d' 

91 91 



X^-y lA + 



1532;^ - 383a;2 + 245x^ - 27x^ 



18(1 



+ 



9r9r rq 
d -d' Jii ' 

9t9t 



-16 + 64a;, 



36a;2 - 93a;3 + 84a;^ 



9(1-^.) 



— - h\x„ 



(9) 
(10) 



^i9t.9tni,,+gi-gim,)F^^'\x^) + 2 (x J , 
'2{9t9(rn, + gi4m,,)F^'^{x^) + 2g(gim^Fo{x^) , 



(11) 



° ^""^ ^ 24(1 - x)3 4(1 -x)4 

(12) 

^ , , -20 + 31x - 5x^ 2x - 3^2 , 

= 12(1-.)^ - 2(13^ • 

This amplitude also contains terms which are divergent, but its magnetic part is finite, as in the 
V — Z case. The fl term has been evaluated in unitary gauge before [2-4] , and our result confirms 
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those given in Refs. [2, 3]/ but differs from that in Ref. [4]. The magnetic, F^-^, terms have also been 
evaluated in unitary gauge previously in Refs. [4, 5], and the resulting expressions are consistent 
with ours, at linear order in g^. 

For V being a gluon, g^^, the diagram in Fig. 1(b) is absent, and we get 



where 



9l 



9s 



167r2 



d' 



m 



w 



{er-^k>^) {qIPl + 9kPK) + ^^'"'KiGlP^ + GIP^) tj , (13) 



gjgj q _ gj gt 

9i9i^^ 4 



l^x^-llxl-xl 8 - 32,Ty + 18x2 
3(1- -J' 3(1- -J' 



(14) 



Gl = 2{gi-g<[m^,-rgi-gim^)Gf'{x^) + 2 Go (^ J , 

= 2{gi-g'[m, + gi-gim,)Gf^{x;) + 2^f (^r^, Go(a;J , 



(15) 



with 



Gfix) 



-2x - 5a;2 + x^ 
%{l-xf 



3a;2 Inx 
A{l-xY ' 



GJx) 



—4 — X — X 



?>x lux 



4(1 -x)2 2{l-xf 



(16) 



Unlike in the V — Z and 7 cases, this amplitude has no divergence. Our expressions for the 
terms containing gl and G^ in the SM limit agree with the corresponding unitary-gauge results 
of Ref. [6], except for the relative sign between them. The magnetic part containing Gq has been 
previously calculated in unitary gauge in Ref. [7], but the result therein misses a factor of |. 

Before proceeding to evaluate the box diagrams, we note that the divergent terms in A^^^^j/^* 
and M.'^^^^i^* above depend on x^. This implies that in the SM limit. 



9t 



V2 



9 

V2 



9i 



9^ ^ 0, 



(17) 



these divergences cannot be eliminated by the GIM mechanism. We have checked that such di- 
vergences would still be present in the absence of the second diagram in Fig. 1(b), which has a 
greater degree of divergence than the others. In contrast, the corresponding terms calculated in 
gauges are finite [8-10]. The k^j'^— J(:k^ parts of the V = 7 amplitudes are, therefore, gauge 
dependent. For physical processes, such as dd' — > uP and dd' — > i^i", we will show later that 
the Xg-dependent divergences contributed by the d — > d'V* diagrams in unitary gauge are canceled 
exactly by the x^-dependent divergences in the box-diagram contributions. The magnetic (a'^") 
terms in the d d'V* amplitude are, on the other hand, gauge independent. In particular, our 
expressions for the and Gl^ terms (at first order in g^) agree with those derived using 
gauges [9-11] 

We also note that the divergent terms in -Mli^diz* ^^^^ ^^^^ present in the amplitude for the 
physical decay Z — > dd'. In that case, /c^-dependent terms need to be kept and, as a consequence. 



^ A factor of Xf seems to be missing from the T*^'^ term in the Eq. (25) of Ref. [3] 
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there are additional divergent contributions. The divergent part of the d — > d'Z* amphtude in the 
/c^ 7^ case is then 



M 



g,(div) 
d-^d'Z* 



9^ ;jfj*\ aUt 



clk^ 



m 



w 



(I6rr4 + fc^)c^fc^ 
3m^ 



+ 



9r9r 



elk' , \ {I6ml + k')cie 
+ 7 Lc. - — 



m 



w 



3m^ 



(18) 



where the fc^-dcpcndcnt terms have been supphed by the second diagram in Fig. 1. Accordingly, 
setting k'^ — rriz — ^w/^w^ finds for Z — > dd' 



M 



g,div_ 
Z-^dd' 



9^ 



1927r2 ci 



d'A-9i gt (164 + 1) + 9i gi (244 X, - 164 - 1) ^rI d. (19) 



Clearly, in the SM limit the m^-dcpcndcnt part of the divergence in the Z dd' amplitude vanishes, 
but the remaining x^-independent part is removed only after applying the GIM mechanism. This 
is in accord with what is found in the literature [12]. 

In evaluating the box diagrams, we assume that the W has SM couplings to leptons and set the 
masses and momenta of the external fermions to zero. Thus, from Fig. 2(a) we obtain 



q,box 



g"^ gj gt 



+ 



1 d -d' 

g gjgj ^ 

1287r2 



{x^ + - 6) A + (x,, x^) dYPRd PTm^l^ , 



(20) 



■^^dd'^e+e- 



g^ gt gj 

1287r2 mly 

1 d -d' 

g gigK 



+ 



1287r2 



(6-xjA--Bi(x,,0) 
(6-xjA-S2(a;^,0) 



dYPLd£%PJ 
dYPnd^PJ , 



where 



3 3{x + y-xy) {Ax'' - 8x^ + x"") In x {Ay' - Sy^ + y^) \ny 

2^"" (l-x)(l-y)+ (y-x){l-xr + (x-y){l-yr 

9{x + y — xy) (4 — x)'x' Inx (4 — y)'y' In y 



+ 



(1 - x)il -y) (y- x)il - x)' (x - ?/)(l - y)' ' 



(21) 



(22) 



and i in Eq. (20) refers to the charged lepton in the loop. Since the divergences in these box 
contributions depend on m^, they cannot be removed by the GIM mechanism in the SM limit, 
Eq. (17). This is different from the corresponding contributions computed in the gauge, which 
are finite [9] 



6 



Combining the dd! Z* ^ vi/ amplitude derived from Eq. (4) and the box-diagram amphtude 
in Eq. (20), we then find 



dd'—*vv 



+ 



327r2 



m 



w 



2 r! -rl' 

327r2 



x, — 6 



2^.+ 4 



dYPRdP^^P^u, (23) 



where 



x(x + 2) 3x(x-2) 
= 8(^+8(^^^""' 



X(x) = 2x- Inx- 4X^{x) . (24) 



X 



In writing down this amphtude, we have neglected the contribution of the Z-magnetic dipole terms, 
with Zl ^, in Eq. (4), which are small compared to the zf^ ^ terms due to the m^y suppression. We can 
see that in the SM hmit, Eq. (17), the divergent part of M.'^^,^^- no longer depends on because 
of cancelation between the m^-dependent divergent terms in the Z-penguin and box contributions. 
The remaining divergence will then disappear when the GIM mechanism operates. 

For dd' £~^£~, the amplitude receives not only the Z- mediated and box contributions, but 
also the photon-mediated one derived from Eq. (8). As a result 



g-gi-gt fSA _ ^^^^^^^ y^^^P^dh^P^L 



327r^ 



+ 



+ 



+ 



+ 



327r2 mfy 

2 d -rl' 2 r 

9 9r9r '< ^ 
i9^ sj 



8A 



+ 8Z,{xUdYPLdi% 



--2x\A + Y{x^) 



dYPRdi^PJ 
dYPRdh/ 



4x,-- A + Z(xJ 



d'a>^^K{FlP^ + FlP^)dh/, 



(25) 



where 



n(^) 



+ 



?>x' 



x{x — 4) 

8(a; - 1) ' 8(x-l)2 



Inx , 



18x^ - 163x=^ + 259^2 - 108x 24x^ - Gx^ - 63^2 + 50x - 8 , 
Zq{x) — ^-TTT TTZ 1 ttt. Inx 



(26) 



144(x - 1)3 



72(x- 1)^ 



5x — 2x^ 

Y{x) = -2x + lnx + 4Yf^{x), Z{x) = 2x - 4x Inx + 8Zf^{x) , (27) 

and we have again neglected the contribution of the Z- magnetic dipole terms in Eq. (4). In the 
SM limit, the divergent terms in the total amplitude are also independent of m due to cancelation 
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among the m^-dependent divergent parts of the penguin and box contributions. This, along with 
the similar cancelation in the dd! — > uV case, explicitly confirms the result of unitary-gauge analysis 
in Ref. [13]. 

In terms of the couplings k^'^ parameterizing the anomalous interactions of quarks and the W 
boson considered in Ref. [1], we have 



9t - ^^.(i+'^y , 



(28) 



where the k's are assumed to be small compared to unity. It is then straightforward to arrive at 
the effective Hamiltonians for dd! — vv and ddi! within the SM, as well as those induced 

by the anomalous couplings, given in Ref. [1]. 

For dd^ dd' , we derive the amplitude from the two box-diagrams in Fig. 2 with quarks d 
and d' ^ d m. the external legs and quarks q and q' in the loops. In this case we adopt the 
parametrization in Eq. (28). With the contribution of the anomalous couplings included to second 
order in it follows that 



•'^dd'-^dd' 



rA \L \L 

2567r2 m^^^ 



+ 



9^ ^q^q' 



2567r^ 



(6 - - Xq,)A - Xq) 
(6 -Xg -a;^,)A - B^{xq,x^) 



d[rPLdld',^aPLd2 



X d[rPLdi d'^lAd2 + d\rPKd, d'^l^Pi^d^) 



9 W'mqm^, 
647r^ 



-|- B^(^Xg, Xg,)^ (^'^^d'^q'd diP^di d2P^d2 + iifd''^^*d' d'lPi^di d^Pi^d^ 



+ {d[ ^ 4) , 
where we have distinguished the two (i*^'^'s, 

= Ajl + «:^;)(l + «:y , 



\ — \'^qd''^qd 1 



(29) 



(30) 



xy — X — y — 2 (4x - 2x^ + x^) In x [4y - 2y^ y^) In y 



{\-x){\-y) {y-x){\-xy 



{x-y){\-yf 



(31) 



and we have neglected finite terms quadratic in containing d! P-^^^dd! o ^^P-^^^d whose con- 
tributions to the neutral-meson mixing would vanish in the vacuum saturation approximation. It 
is also straightforward to obtain the corresponding effective Hamiltonians given in Ref. [1]. This 
dd! — > dd amplitude in unitary gauge has been previously calculated in Ref. [14], but our result for 
the terms disagrees with theirs. 

Finally, we remark that some of the physical amplitudes we have found above generated by the 
anomalous couplings contain divergent parts. As discussed in Ref. [1], this is due to the fact that 
the effective theory with anomalous couplings in Eq. (1) is not renormalizable, and the divergences 
are understood in the context of effective field theories as contributions to the coefficients of higher- 
dimension operators, which are not included in our analysis. 
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In conclusion, we have provided the detailed derivation of loop formulas summarized in our 
recent work evaluating several two- and four-fermion FCNC transitions induced at one-loop level 
by the anomalous charm-H^ couplings. In treating the loop diagrams, we have worked in unitary 
gauge and used dimensional regularization to handle the divergences arising in some of the loops. 
We have compared our results with the corresponding ones obtained in the literature using unitary 
and gauges, where available. In particular, we have discussed subtleties in the cancelation of 
divergences from individual contributions to some of the physical amplitudes and derived expressions 
consistent with those obtained in gauges. Finally, we have provided, perhaps for the first time, 
the unitary-gauge expressions for the separate diagrams contributing to the d d'V* amplitude. 
These particular results can be applied to other models involving fermions and gauge bosons with 
interactions similar in form to those considered here. For example, the contribution of the Z' boson 
to the b ^ transition can be easily found from the magnetic terms in our general expression for 
the d — > d'V* amplitude, with appropriate couplings. 
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APPENDIX A: AMPLITUDES FOR d d'V 



The Lagrangian describing a neutral gauge-boson V coupling to a fermion-antifermion pair, //, 
or a ly-boson pair can be expressed as 

£v = -h'{L)P^ + R]Pn)f^, + wl[V'{w;,W-^ - W+^W-J) - V^^W^W-] , (Al) 

where Lj, and g^^ contain the coupling constants, and X^i, — dfj,Xi, — duX^. This and Cudw 
in Eq. (1) provide the relevant vertices for the d d'V* transition at one-loop level, leading to 
the diagrams in Fig. 1, where the external quarks d and d' are on-shell and V is off-shell. The 
resulting amplitudes for the contributions (a,b,c), with q being the quark in the loops and V having 
a four-momentum k, are, respectively. 



w{a)<? 



d -d' 

167r2 



'rY-^ly 



■ X, 



fv ^ jl^ S ' fv 



m 



w 



+ V, (x J + R^^ V, (x J ) (m,, Pl + m, P^) 

4^ d 

+ {Ql ^ 9r, L^, ^ R^, Pl ^ Pr) , 



m 



w 



(A2) 
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ia^'^e* k 



i(j>^''e* k 

+ ^ioW^,(^l^r'^l + ^l ^r^r)— , (A3) 



= - -^2w)^r(5^:5^:'^^^L+^R^R'^x)^ , (A4) 

where A is given in Eq. (3), e and A; are the polarization and four- momentum of V, respectively, 



-gi - {giy , -gi = {g^y , ^/ = ^ , <^"' = , (A5) 



N 3 ^ , , -7a; + a;^ 4a; - + , 

V,{x) = ---x + 2xlnx, = 2(1 -x)^ 

, 2 10a; - 5a;2 - lla;3 (2 - Sa;)^ , 

Voix) = 1 ; — — mx , 

' 9 36(1 -a;)3 6(1 - a;)^ 

, , 1 llx - 7a;^ + 2a;^ 2x \nx 

^4(^) = o+ KT. -T-, + 



3 9(l-a;)3 3(1 - a;)^ ' 

= 3 + 12(1 + 5031)1 . = 4-4 ''^<^' 

^ ^ , , 9a; — 3a;^ 3a; In x 

Vj{x) = -1 



4(1 -a;)2 2(1 -a;)3 ' 
, 25 27a;-113x2 + 83a;3-9a;^ 31^3 - 28x^ + 3a;5 , 

Vq(X) — 1 ; — ; lu a; , 

^ 18 24(1 -a;)3 12(1 - a;)^ 

OJU ^JU OJU 111 Jj ^ ^ / ■. 

^'^""^ " 6 " 4(i_^)3 + 2(i_^)4 ' ^io(^) = 6y4(a;) + ^^(a;) + 3V,{x) - 2V,{x) . 

In obtaining these results, wc have assumed that k and rrid^d' are small compared to m^^, kept terms 
to second order in k, and taken at the end of the calculation the usual limit of vanishing rrid^d' for 
the nonleading terms [8, 9]. It is worth noting that the x^-dependent part of the divergent term in 
each of these amplitudes arises from the p^p^ terms of the W propagators. 

The separate expressions above for M-^^-^^v, in unitary gauge may have been written down for 
the first time in this work. They are general and follow from any effective interactions of the form 
given in Eqs. (1) and (Al), provided that the masses of the external fermions and the momentum 
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of the external gauge boson are much less than the gauge-boson mass in the loops. In the case of 
the standard model, the couplings of V — Z, 7, g^^ in Cy are parametrized by 

(AT) 



L] ^ R] ^ eQf , Lf = Rf = g,t, , 



9w ^ , 9w = ^ = 9s^ , = ' 

= cos0^ , = sin6lw , 

where I^f and Qj are as usual the third component of weak isospin of / and its electric charge, 
respectively, and i„ are the color-SU(3) generators satisfying Tr{tatb) — |(^a6- 
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